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1. Introduction 



In the study of medium to heavy mass nuclei, nuclear Density Func- 
tional Theory (DFT), based on the self-consistent Hartree-Fock-Bogoliubov 
(HFB) method and its extensions, is the theoretical tool of choice. As new 
exotic beam facilities are being built or proposed to be built around the 
world, DFT is on the edge of becoming a predictive theory for all nuclei 
but the lightest. This is not only true for ground state properties, such as 
binding energies, radii or multipoles of the density, but also for low-energy 
spectroscopy and decay probabilities. 

These advances are possible thanks to the development of better energy 
functional and to the increase of computer resources. However, the needed 
accuracy and predictive power for unknown regions of the nuclear chart 
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still loave a lot of room for improvements. Thus, DFT faces important 
challenges in its quest for the truly universal energy density functional 
which, in addition to allowing reliable calculations of extended nucleonic 
matter, nuclear masses, beta decay rates or fission lifetimes, is crucial to 
the understanding of tlio astrophysical r-process relevant to the nucleo- 
synthesis of heavy elements. 

In the present contribution, some of those challenges are discussed**. As 
a matter of illustration, preliminary results of an ongoing study address- 
ing the idea of using new ab-initio results to adjust the Skyrme energy 
functional are presented. 

2. Challenges 

The last decade has seen tremendous achievements in developing micro- 
scopic methods to describe medium to heavy nuclei. The most striking ones 
have come from nuclear DFT based on the self-consistent HFB method and 
its extensions (for a recent review, see Ref.^). Such methods rely on reahstic 
effective interactions/functionals (either non-relativistic Skyrme or Gogny 
functionals, or relativistic Lagrangians) as the only (phenomenological) in- 
put. Specific correlations beyond the static mean field are now (almost) rou- 
tinely incorporated through restoration of broken symmetries and configu- 
ration mixing calculations (projected-Generator-Coordinate-Method DFT). 

Beyond the successes encountered in the development of the method 
and in the comparison with experimental data, several formal and practical 
challenges still have to be overcome in order to obtain a unified descrip- 
tion of all systems throughout the chart of nuclei, in particular with the 
perspective of the construction of new radioactive beam facilities like RIA. 

2.1. Long term strategy 

Energy functional practitioners need a better-defined long term strategy in 
the spirit of the "Jacob's ladder" which is referred to in electronic DFT.^ 
Indeed, constructing and adjusting functionals is a tricky game where it 
has been difScult to learn from other groups' work. Each group uses its 



''This is not an exhaustive list of challenges. In particular, it does not cover all meth- 
ods related to DFT, i.e. QRPA. Some elements discussed presently reflect the overall 
view of the DFT community and constitute the outcome of the workshop "Towards a 
universal density functional for nuclei", September 26-30, 2005, Institute of Nuclear The- 
ory, Seattle, USA. Of course, some ideas discussed here reflect our personal view on the 
matter. 
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own construction/fitting strategy, itself evolving in time in a way which is 
sometimes not systematic enough to allow for significant progress. 

First, functionals at different levels of the many-body treatment must be 
consistent. For example, the HFB and the projected-GCM functionals must 
be connected in a clear way, the latter including the former as a particular 
case. Thus, projected-GCM DFT should be motivated from first principle 
and the HFB level derived from it. This is still lacking. Second, any im- 
proved/complexified functionals should include all features validated by the 
previous/simpler ones. Third, a list of (truly) benchmark problems should 
be argued upon and used to make the energy functional consistent with 
ab-initio calculations of systems which are essentially solvable exactly; i.e. 
Brueckner-Hartree-Fock (BHF)^ and variational chain summation (VCS) 
methods"' in infinite matter, Green Function Monte Carlo (GFMC)^ method 
and No Core Shell Model (NCSM)*' for light nuclei, later followed by Cou- 
pled Cluster (CC) methods''' for heavier systems. An example of that wih 
be touched upon in section 3. Fourth, practitioners need to make progress 
regarding fitting algorithms. Systematic errors must be estimated and co- 
variance analysis of parameters need to be performed to determine their 
relevant combinations.^ Also, one has to assess the improvements brought 
about by further ingredients: new theory input and/or additional experi- 
mental data, i.e. superdeformed states, fission isomers/barriers of (exotic) 
nuclei. Finally, one can test new fitting algorithms (re-annealing,^ genetics) 
to make sure that the absolute minimum in parameter space is reached 
during the fit. It is only fohowing such a "strict" strategy that significant 
progress will be made in the near future. 

2.2. Improved phenomenology 

Improving single-particle energies coming out of the functional might be 
one of the most important subject of focus for the near future. Indeed, 
many unsatisfactory single-particle energy spacings provided at the HFB 
level spoil features of low-energy spectroscopy which cannot be corrected for 
by including further correlations through beyond mean-field techniques.^'' 
This has not only to do with the shell evolution towards the drip-lines but 
also with single-particle energies in stable systems. In any case, it is clear 
that data on particle /hole states in very neutron rich nuclei will help by pro- 
viding shell evolution along long isotopic chains. For instance, information 
on particle/hole states around "^^Ni should be accessible with RIA. 

However, there are issues as far as the fitting strategy is concerned. 
Indeed, one cannot fit single-particle energies directly to experiment with- 
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out including their renormalization through the coupling to surface vibra- 
tions.^^ One would need to adjust functionals through beyond mean-field 
calculations of odd nuclei, which is too time consuming to be done in a 
systematic way within projected-GCM DFT. One could use other methods 
to estimate the effect. 

One has indications that the previous problem might be partly related 
to the omission of the tensor force in DFT. The latter could explain a lot 
of closing/opening of sub-shells as one fills major shells. This is related to 
the attraction/repulsion between neutron (proton) j = I + 1/2 and proton 
(neutron) j — I — 1/2/1 + 1/2 orbits provided by the tensor interaction. 
The = 32 sub-shell closure in neutron rich isotopes (^^Ti) is one possible 
example. Several groups are now working on implementing the tensor force 
in the Gogny and the Skyrme functionals. ^^"^^ It is important to note that 
the tensor force has been mostly disregarded in DFT so far because of a lack 
of clear-cut experimental data to adjust it. Again, observation of (relative) 
single-particle energies in very exotic nuclei could be of great help in that 
respect. 

Masses, separation energies, densities, deformation, individual excita- 
tion spectra and collective excitation modes such as rotation or vibration, 
depend significantly on the superfiuid nature of nuclei. The role of pairing 
correlations is particularly emphasized at the neutron drip-line where the 
scattering of virtual pairs into the continuum gives rise to a variety of new 
phenomena in ground and excited states of nuclei. 

Despite its major role, our knowledge of the pairing force and of the 
nature of pairing correlations in nuclei is rather poor. The structure of the 
functional as far as gradient versus (isovector) density dependencies are con- 
cerned has to be clarified. ^'^"^^ Proper renormalisation techniques of local 
functionals have been recently proposed but their content must be further 
studied. Also, particularly puzzling is the situation regarding beyond- 
mean-field effects associated with density, spin, isospin fiuctuations.^-'^'^^ 
Finally, the infiuence of particle-number projection and pairing vibrations 
has to be characterized through systematic calculations. 

At this point, all modern phenomenological functionals provide similar 
pairing properties in known nuclei, i.e. odd-even mass staggering (OES) 
between ^^^Sn and ^^^5n, while their predictions diverge in neutron rich 
systems, where having masses up to ^'^^^^^'^Sn or ^-'^A''^ with an accuracy 
SE = 50 kcV would already allow one to discriminate between most models. 
Pairing can be probed in other ways and data on high spin states, low- 
energy vibrations, two-nucleon transfer cross sections (feasible with beams 
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of 10^ part/s) in exotic nuclei will be complementary and all very valuable. 

2.3. Connection to underlying methods 

The energy functional has to be rooted into a more fundamental level, even- 
tually connected to QCD. Doing so is not only important from a heuristic 
point of view but also from a practical one. Indeed, the main problem of 
current phenomenological functionals is the spreading in their predictions 
away from known nuclei. At tho samo time, phenomenology indicates that 
gradient and/or density dependencies of usual Skyrme functionals are too 
schematic/limited (see section 3.2). Thus, we need guidance beyond a fit 
on existing data. Hopes are put into deriving functionals from Effective 
Field Theory (EFT) which has been successfully applied to few-nucleon 
systems and which offers a comprehensive picture of the relative impor- 
tance of three-body and two-body forces. This could be done with the 
use of the renormalization group method in the spirit of Viowk^^ from which 
a perturbative many-body problem seems to emerge in infinite matter.'^'' 
The connection with finite nuclei will call for methods such as the Density 
Matrix Expansion. 

2.4. Grounding beyond- mean- field DFT 

Ultimately, the functional has to incorporate, in a microscopic manner, long 
range correlations associated with large amplitude vibrational motion and 
with the restoration of broken symmetries. This is of particular importance 
to finite self-bound systems. The projected-GCM extension of DFT has 
proven to be one of the most promising ways of doing so. For the future, the 
first challenge will be to treat more collective degrees of freedom within the 
GCM, such as the pairing field or the axial mass octupole moment. So far, 
only the axial quadrupole moment has been considered in systematic GCM 
calculations.^ Second, configuration mixing and symmetry restoration will 
have to be properly formulated, from first principle, within the context of 
DFT. While one can use guidance from the case of the average value of a 
Hamiltonian in a projected state, DFT presents key differences which have 
to be considered seriously. Because of the lack of theoretical foundation, 
one has rehed on ad-hoc prescriptions so far to extend the HFB functional 
to the projected-GCM level, in particular as far as density dependencies 
in the functional are concerned. Even worse, particle-number projected 
DFT is ill-defined as was recently reaHzed^^ and fixed. ^° 
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3. The Skyrme energy functional : use of ab-initio inputs 

In relation to section 2.1, we now briefly present a recent attempt to re- 
flne the SLyx family of Skyrme functionals^^ by using additional ab-initio 
inputs from calculations of infinite nuclear matter. Our goal here is to high- 
light the general lessons one can learn from such an attempt. 

3.1. Constraining the isovector effective mass m* 

Our original goal was to study the effect of the neutron-proton effective 
mass splitting m* — m* with isospin asymmetry I = {pn — Pp)/p on the 
properties of the Skyrme functional. It is of importance since the effective 
mass controls the density of single-particle energies at the Fermi surface, 
and hence influences masses and shell correction in superheavy nuclei as 
well as how static and dynamic correlations develop. 

While the isoscalar effective mass m* (~ 0.8) is well determined by 
the ISGQR,"^'^''^"^ constraints on the isovector one m* (w 0.7 — 0.9) via the 
IVGDR are not decisive enough. On the other hand, recent BHF calcula- 
tions of asymmetric inflnite nuclear matter, with/without 3N force, predict 
Aml_p w 0.22 for / = 1^4 ( Am*_p « 0.13 for Dirac BHF^b.sb). 

Although the amplitude of the splitting cannot be considered as a ab- 
initio benchmark, its sign is more deflnite. Thus, we decided to construct 
Skyrme functionals with different values of Am*_p with the idea of im- 
proving their isovector properties by using predictions from ab-initio calcu- 
lations. From that point of view, it is interesting to note that the SLyx pa- 
rameterizations were fltted to the Equation of State (EOS) of Pure Neutron 
Matter (PNM) with the idea of improving isospin properties of the function- 
als. One consequence was to generate functionals with Am*_p < 0, in op- 
position to ab-initio predictions. On the other hand, older functionals such 
as SIIP"^ and SkM*,^^ which were not fitted to PNM, had Am^.p > 0. 
The same exact situation happens for the Gogny force. Thus, improving 
global isovector properties (EOS) seems to deteriorate state-dependent ones 
(rn*) with currently used functionals. 

In any case, we tried to adjust three parameterizations (/_, /o, /+) 
using the same fitting protocol as SLyb but with different values of m*. ^. 
Their bulk properties in SNM are summarized in the table below. They are 
all identical and very close to SLy5', except for the effective mass splitting. 



'^The only difference in the protocol is that static spin-isospin instabilities below p = 
2psat are excluded through the use of Landau parameters. SLy5' is as SLy5 except for 
the ab-initio EOS used in the fit; see Fig. f . 
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Psat 


E/Asat 




a-i 


* 


Am*_j, 


SkM* 


0.160 


-15.770 


217 


30 


0.79 


0.356 


SkP 


0.162 


-15.948 


201 


30 


1.00 


0.399 


SLy5' 


0.161 


-15.987 


230 


32 


0.70 


-0.182 


/- 


0.162 


-16.029 


230 


32 


0.70 


-0.284 


fo 


0.162 


-16.035 


230 


32 


0.70 


0.001 


f+ 


0.162 


-16.036 


230 


32 


0.70 


0.170 



First, a few lessons were learnt. Raising m*/m to 0.8 is difScult because 
of spin-isospin instabilities. With ml/m, = 0.7, it is impossible to reproduce 
the high density part of PNM EOS and fix Am* _p > at the same time'^. 
It explains why the SLyx forces predict the wrong sign of the effective mass 
splitting in neutron rich matter. Using two density-dependent terms in the 
functional cx pj^'^; Po^^"'^ allowed the construction of (/_, /o, /+). 

In addition, 5Ly5-like functionals with Am*_p > generate new finite- 
size isospin instabilities which become more dangerous as Am* _p becomes 
more positive. The instability tends to split neutron and proton densities by 

gaining energy through the term ac Cf''' [v pi^ in the functional, where 
pi is the scalar-isovector density. The table below gives typical values of 
Ci''. We could check that SkP^^ already presents such an instability and 
that, empirically, must be greater than Ri — SOMeF to have an insta- 
bility free functional. Those finite size instabilities need to be quantified 
and better controlled in the fit of the functional. For such a purpose, we 
have extended the control to finite wave-length instabilities by calculating 
the RPA response function in symmetric nuclear matter (SNM).^°'^^ 





f- 


SLy5' 


fo 


/+ 


SkP 


Am*_p 


-0.284 


-0.182 


0.001 


0.170 


0.399 




-5.4 


-16.7 


-21.4 


-29.4 


-35.0 



The previous discussion shows already the type of problems /information 
arising from our attempt to improve on the fitting protocol of the Lyon 
functionals by using more ab-initio benchmark inputs. Now, the left panel 
of Fig. 1 shows PNM and SNM EOS as obtained for (/_ , fo, /+, SLy5') and 
as predicted by variational chain summation methods." The asymmetry 
energy a/ as a function of density is displayed on the right panel. The 



'^For all formula of relevant quantities, we refer the reader to.' 
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message to extract from such a comparison is that the four forces have 
the ability to reproduce both microscopic EOS with the same accuracy. 
However, is that enough to characterize the spin-isospin content of the 
functionals and assess their quality? In order to answer such a question, we 
use further ab-initio predictions and look at the contributions to the SNM 
EOS per spin-isospin {S,T) channels. 




0.08 0.16 0.24 0.32 0.4 0.08 0.16 0.24 0.32 0.4 

P [fm"^] p [fm-^] 

Figure 1. Left : SNM and PNM OES as given by functionals discussed here (see text), 
compared with VCS results by Akmal et at.^ Right : isospin symmetry energy vs. density 
for the same set of functionals, compared to variational calculations by Lagaris and 
Pandharipande.** 



3.2. Separation of the EOS into (S,T) channels 

In this section, we discuss the contributions to the potential energy'^ of SNM 
from the four two-body (S, T) channels. We compare our results with those 
predicted by BHF calculations using Argonnc V18 two-body interaction and 
a three-body force constructed from meson exchange theory. By virtue 
of antisymmetrization, each spin-isospin channel is linked to a given parity 
of the spatial part of the effective force. However, and as motivated by the 
DME, the Skyrme functional is to be seen as the average of an effective 
vertex over the Fermi sea.^^ In particular, the angular dependence of the 
original nucleon-nucleon force have been averaged out. Thus, the fact that 
only S and P waves arc explicitly considered in the Skyrme functional is 
illusory. In any case, one has access to the averaged spin-isospin content by 
looking at the different (5, T) contributions to the energy. 



'^The difference between the total energy and the kinetic energy of the corresponding 
free Fermi gas. 



4th February 2008 9:33 



WSPC - Proceedings Trim Size: 9in x 6in riacproc 



Using projectors on spin/isospin-singlet and spin/isospin-triplet states, 
it is possible to split the potential energy Ep into {S, T) contributions E^'^ . 
In SNM, one finds : 



i: (0.(1 + x„.)/«'= + j^t,(l + X,) C^)'" p'". (2) 



160 ' ' V 2 
16 

El ^3^^ . 9 . ./3.2\^/^ 

^ii = ^..fl..oir^V^^/3 



where {ti,Xi) are usual coefficients of the Skyrme functional, whereas 
(toi, xoi) characterize the density-independent zero-range term and the two 
density-dependent ones. 



> 



3 
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o -10 

fi -15 
-20 
5 -25 
U 













SLy5' 



5 ^ 

-10 S3 

-5 > 

-10 § 

-15 ^ 

-20 3 

-25 H 

-30 < 



0.1 0.2 0.3 



0.1 0.2 0.3 0.4 



53 



Figure 2. Energy per particle in eacii {S,T) channel for SNM, as a function of density. 
Crosses refer to the BHF calculations.*' 



Results are plotted against BHF predictions in Fig. 2. First, one can 
observe that the results are rather scattered. Second, the main source of 
binding, from the {S,T) = (0, 1) and (1,0) channels, is not weh described 
and the saturation mechanism is not captured. It is thus clear that, even 
though all four forces reproduce perfectly the PNM and SNM EOS, they 
do not have the same spin-isospin content, and that the latter is in general 
rather poor. Thus, fitting the global EOS is an important element but it 
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docs not mean that the functional has good spin-isospin properties. One 
needs to do more and fitting ab-initio predictions for the contribution per 
(5, T) channel seems to be a good idea. However, one needs to make sure 
that the data used are truly benchmarks calculations. This calls for predic- 
tions from other ab-initio methods. 

The most obvious discrepancy appears in (0, 0) and (1, 1) channels where 
Skyrme and BHF data have opposite signs above saturation density. The 
SLy5 parameter set shows a particular behavior in the (1,1) channel due 
to the choice of X2 = — 1 to prevent ferromagnetic instabilities in PNM. 
Note that these two channels are taken care of in the Skyrme functional 
by the density-independent P-wave term only. The upper-right panel of 
Fig. 2 points out the tendency of Skyrme parameterizations to be attractive 
in polarized PNM and cause a collapse of its equation of state at high 
density. At lower densities, the BHF data show a distinctive behavior, being 
slightly attractive below po and repulsive above. This feature cannot be 
matched by the standard Skyrme functional since it exhibits in this channel 
a monotonous behavior with density, whatever the value of {t2,X2)- 

It is also worth noticing that the failure in the (1,1) channel becomes 
more and more prominent as one makes the splitting of effective masses 
closer to the ab-initio predictions (force f4). The effective masses being 
governed by the momentum dependent terms of the force, it is not a surprise 
that the modification of the former impacts the (0,0) and (1, 1) channels. 
What change in the coefficients entering Eqs. (1-4) stems only from the 
variation of m* and the associated rearrangement of parameters in the 
functional, most notably the C^f coefficients closely related to the surface 
and surface-symmetry energies. The quite tight requirements on the latter 
imply that the four parameters of the non-local terms in the standard 
Skyrme energy functional would be dramatically overconstrained if we were 
to add the (S, T)-channel decomposition in the fitting data. 

The rather poor properties of the functional in the (0, 0) and (1, 1) chan- 
nels, the degradation of the latter as the effective mass splitting is improved, 
the finite-size instability problem, the idea of using the ab-initio {S, T) con- 
tributions in the fit, call, at least, for a refinement of the odd-L term in the 
sense either of a density dependence or a higher-order derivative term. The 
latter being prone to numerical instabilities and interpretation problems, a 
density-dependent k' • k term remains as one of the next potential enhance- 
ments to be brought to the Skyrme energy functional. Phcnomenological 
constraints on gradient terms are mainly related to the surface of nuclei, 
i.e. low-density regions. One can expect that, to first order, BHF data in 
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{S,T) = (1,1) channel can be matched with an extended functional while 
retaining a good agreement with other (experimental) data. It is less clear in 
the (0, 0) channel but further exploration of the extended parameter space 
may bring Skyrme and BHF data in better agreement. 

4. Conclusions 

Density Functional Theory faces important challenges in its quest for the 
universal energy density functional applicable over the entire nuclear chart 
with a true predictive power. In this contribution, we discuss some of those 
challenges. Then, we exemplify one of them by showing the use which can 
be made of further constraints from ab-initio calculations and how usual 
Skyrme functionals are too schematic/limited as well as over constrained 
as soon as one tries to bring more microscopic character to them. 
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